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YIELD  AND  BLAST  ANALYSES  WITH  A UNIFIED  THEORY  OF  EXPLOSIONS  (UTE) 


\1 


Summary 


Yield  is  the  most  significant  single  number  measured  on  any  explosion  because 
all  effects  —sympathetic  reactions,  fragments,  blast  damage—  derive  from  it. 

UTE  offers  the  only  adequate  way  to  relate  all  explosions:  non-ideal,  any  media, 
over  all  ranges  in  air,  underwater,  underground,  confined  spaces,  heavy  cases  etc. 
The  Form  Factor  and  Lead  Time  are  new  extensions  intended  for  safety  analyses. 

A key  idea  in  the  form  factor  is  to  define  average  energy  density  in  the  blast  wave 
relative  to  the  peak  value  at  the  shock;  it  is  the  tacit  assumption  in  scaling  now. 
Lead  Time  means  the  difference  in  TOA  between  a sound  signal  and  a shock  wave; 
it  scales,  is  a sensitive  measure  of  yield  and  is  nearly  constant  at  long  range. 
Applications  include:  absolute  measure  of  prompt  and  delayed  yields  for  blast, 
sympathetic  reactions,  fragments,  propellant  yield,  surface  effects,  analysis 
with  sparse  data  and  simplified  instrumentation. 
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YIELD  AND  BLAST  ANALYSES  WITH  A UNIFIED  THEORY  OF  EXPLOSIONS 

Francis  B.  Porzei 
Naval  Surface  Weapons  Center 

1.  UTE  Methods  for  New  NESIP  Problems 

Yield  means  the  blast  or  hydrodynamic  energy  released  by  any  explosive 

and  is  the  single  most  significant  number  to  be  measured  in  tests  of  any  explosion. 

All  effects  derive  from  yield  and  in  principle  can  be  predicted  using  it: 

* Primary  fragments:  their  sizes,  shape,  number,  initial  velocity  of  trajectory, 

* Close-in  blast,  regarding  both  initial  containment  and  secondary  fragments, 

* Low  pressure  damage,  especially  to  specify  the  range  where  1 psi  occurs, 

* All  sympathetic  reactions  -detonation,  deflagration,  burning,  their  degree- 
be  it  directly  from  blast  or  via  fragment  and  thermal  loads  it  produces. 

In  fact,  the  Maximun  Credible  Event  (MCE)  really  means  the  overall  yield. 

The  Unified  Theory  of  Explosions  (UTE)^  was  developed  for  just  such  needs. 

UTE  offers  the  only  general  way  to  describe  any  explosion:  nuclear,  non-ideal  HE, 
over  all  ranges  in  air,  underwater,  underground,  confined  spaces,  heavy  cases  etc. 

UTE  offers  two  dozen  concepts  as  tools  to  treat  dozens  of  real  non-ideal  effects, 
nearly  all  being  unknown  or  ignored  in  idealized  classic  theory  and  hydrocodes. 

A key  idea  in  UTE  is  "prompt  vs  delayed"  energy;  it  asserts  that  natural  processes 
release  some  energy  instantly,  some  more  slowly,  reinforcing  blast  farther  out; 
much  is  trapped  behind  the  negative  phase  too  late  ever  to  support  the  shock  front. 
This  separation  is  manifest  by  phenomena  like  afterburning,  secondary  shocks, 
and  the  most  dramatic  feature  of  any  explosion:  The  fireball  _is  delayed  energy. 

2 

The  NESIP  Technology  Base  itself  rests  on  the  Unified  Theory  of  Explosions  and 
much  success  in  NESIP  tests  is  due  to  versatile  and  accurate  analyses  with  UTE. 

Current  NESIP  problems  now  raise  new  and  more  specific  questions  about  yield. 

For  sympathetic  reactions  and  in  the  design  of  inhibitors  to  prevent  them: 

* How  much  energy  is  released  in  a partial  or  in  a low-order  detonation? 

* What  are  the  actual  prompt  and  delayed  fractions  in  afterburning  explosives? 

* How  to  live  with  the  large  scatter  of  pressures  measured  in  the  real  world? 

* How  to  live  with  the  narrow  range  of  feasible,  affordable  measurements? 

* What  are  the  absolute  energies  involved  in  various  modes  of  energy  release? 

For  hazards  involving  propellants: 

* What  are  the  yields  of  a propellant  on  an  absolute  basis,  detonated  alone? 

* If  set  off  by  an  explosive  warhead,  how  much  does  the  propellant  add 

to  the  prompt  or  delayed  yield  of  the  main  explosion?  Any  new  hazards? 

To  meet  these  new  needs  for  NESIP,  two  major  advances  have  been  developed  for  UTE: 

* A "form  factor"  method  for  bookkeeping  the  energy  within  a blast  wave, 

to  describe  variable  rates  of  afterburning,  notably  in  heavily  aluminized  HE. 

* A "lead-iime"  method,  a simple  reliable  way  to  get  yield  from  time-of -arrival. 

They  have  widespread  application  to  many  explosion  problems  for  NESIP  and  others. 

The  purpose  of  this  paper  is  to  describe  these  new  methods  briefly,  show  the  code 
and  to  test  them  by  comparisons  with  a broad  spectrum  of  field  measurements. 

1,  All  references  are  listed  at  the  end  of  the  text,  before  Appendix  A. 
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2.  Form  Factor  Concept  and  Method 


The  form  factor  F is  defined  to  mean 
the  average  hydrodynamic  energy  W+K  in  the  wave 
relative  to  their  peak  sum  E at  the  shock  front. 

It  is  the  ratio  of  the  area  sfiown  here  as  F 
to  the  "square  wave"  as  it  E were  constant  at  fc.?. 

If  we  define  the  yield  Y{R)  at  any  shock  radius  R as  the  inteprated  sum 

Y(R)  = 4pi  / (W  + K) r dr 
0 

normalize  the  integral,  multiplying/dividing  by  the  peak  overpiessure  Ps-PQ, 

and  by  the  mean  value  of  any  factor  (like  the  available  fraction  A of  energy  at  P) 

Y (R)  = 4p  i R3(F  -P  )A  / rW  + K-,3,-(r)2dr1 
3“  S 0 / [P-”P]A[(R)  ~R 

0 S 0 

Thus  F becomes  the  dimensionless  fraction  specified  by  the  definite  integral. 
Rigorously,  we  can  just  simply  define  a form  factor  F such  that 
yield  volume  , O'pressure  mean  A form  factor 
Y(R)  = ( 4p  i / 3 ) R3  P A F 

This  definition  for  F is  deceptively  simple  but  is  a powerful  hydrodynamic  tool. 

Many  man-years  and  $millions  were  spent  since  World  War  II  on  elaborate  hydrocodes, 
mostly  to  calculate  pressure-distance  curves  with  highly  over-simplified  models. 

Yet,  both  A and  F are  readily  prescribed  by  the  overpressure  ratio  (P  /P  -1); 
so  we  can  always  obtain  the  shock  radius  R a*.  .;ny  pressure  P simply  rrom 

R3  = Y(R ) P =overpressure,  units  consistent  with  '!  and  R 

( 4p  i 73T  PAF 

To  calculate  R,  the  code  decreases  Y(R)  from  its  initial  .-air'e  Yq  by  decrements 

dY  = -4pi  Q R2dR  + Afterburning 
losses  gains 

This  can  be  done  in  bold  steps,  decreasing  the  pressure  about  25,'o  p*  each  step. 

The  same  steps  are  used  also  for  integrating  the  time  of  arrival  of  the  shock  wave. 

A and  F always  appear  together  and  here  is  an  exact  way  to  calculate  AF. 

For  exposition,  let  G here  mean  the  net  loss  (waste  heat)  and  gain  (afterburning!. 
Then  the  ratio  of  the  dissipation  equation  (dY=)  to  the  definition  for  F (Y=): 
dY=  -4pi  Q R^dR  3(Q/P)dR 
Y = (4p i /3 )RPAF  = AF  R 

gives 

dY/Y  3Q/P  and  AF  = 3Q/P 

6RlR  = AF  din  Y/ din  R 

The  machine  code  uses  the  local  value  of  din  Y/dln  R from  each  previous  step, 
because  din  Y/dln  R varies  slowly,  from  about  -.5  at  high  pressure  to  -1  at  low. 

Thus  AF  is  bounded  between  6G/P  at  high  pressures  and  3G/P  at  low  pressure. 

Figure  1 shows  these  bounds  and  the  transition  region  for  the  function  AF  vs  P. 

2 

Analysis  shows  that  AF  goes  like  A(shock)/3  at  high  pressure,  like  AP  /3  at  low, 
and  a suitable  approximation,  to  a few  percent  in  Y,  without  use  of  din  Y/dln  R, 
and  with  a single  parameter  1/3  for  both  high  pressures  and  low,  is: 

AF  S A( shock ) 

3d  + P0I P)2 
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3.  Lead-Time  Concept  and  Method 


Lead-time  is  a time-of-arrival  method  for  measuring  yield.  With  it  we  seek: 

* An  over-all  measure  of  the  shock  history,  sensitive  tn  its  early  behavior. 

* A way  to  circumvent  uncertainty  in  yield  from  scatter  in  pressure  measurements. 

* Simpler  instrumentation  than  the  sophistication  needed  for  good  pressure  data. 

3 4 3 

Time-of-arrival  is  an  excellent  measure  for  high  pressure  supersonic  blast  ’ ’ ; 

but  at  low  pressure,  TOA  becomes  sonic  and  an  insensitive  measure  of  shock  strength. 
On  the  other  hand,  time-of-arrival  can  be  measured  with  exceedingly  high  accuracy. 
Let  us  then  measure  the  difference  in  time  between  the  shock  and  a sound  signal 
and  define,  at  any  range  R: 

Lead-time  = sound  arrival  - shock  arrival  time 
LT  R/C  - T 

This  quantity  ought  to  and  does  scale,  is  a sensitive  measure  at  low  pressures. 

Best  of  all,  it  becomes  insensitive  to  the  range  at  acoustic  strength,  so  that  one 
does  not  need  an  accurate  gauge  location  —if  the  sound  arrival  is  also  measured. 


As  shown  in  Figure  2a,  the  early  shock  is  highly  supersonic:  U»C  . There 
R = /udt  or  T = /dR/U  are  both  sensitive  measures  of°yield. 

If  we  plot  In  R vs.  In  T as  in  Figure  2b,  sound  speed  is  a straight  line,  slope  1. 

But  the  shock  time-of-arrival  approaches  it,  partly  because  of  the  logarithmic  plot. 
Also,  the  lead-time  ceases  to  grow  as  the  "overvelocity"  vanishes  at  low  pressure. 

As  shown  in  Figure  2c,  the  lead-time  approaches  a scalable  constant  at  long  range. 

In  machine  calulations,  time-of-arrival  adds  up,  using  the  same  steps  as  for  Y(R> 
TOA  = dR/U  U = local  shock  velocity 

Because  of  the  finite  step  size,  an  average  value  for  1/U  is  used;  thus 

"dT"  = 'Ll  -Tj  = . l/U,MRul  • R,V2 

The  time-of-arrival  and  the  lead-time  are  scaled  just  as  for  distance  scaling. 

When  the  shock  is  strong,  it  is  convenient,  and  a more  independent  measure  to  scale 

Relative  yield  = (Measured  TOA/reference  TOA)^. 

When  the  shock  is  weak  (below  the  transition  pressure)  we  scale  lead-times  as 
Relative  yield  = (measured  LT/reference  LT). 


Figures  3 and  4 illustrate  practical  reasons  for  developing  the  lead-time  method. 

Large  scatter  in  field  pressure  measurements  is  no  doubt  real,  probably  intrinsic, 
because  a pressure  gauge  "feels"  only  the  pressure  at  its  surface,  regardless  of 
how  rapidly  pressure  may  vary  in  a boundary  layer  next  to  that  surface®’  . 

Figure  3 shows  how  pressure  suffers  from  real  variations  both  in  space  and  in  time 
Spacewise,  the  dust-laden  boundary  layer,  brush,  rough  terrain  all  deplete  yield. 
Timewise,  the  shock  jets,  and  "rings"  is  it  goes,  in  patterns  that  shift  with  time. 

The  time-of-arrival  will  more  nearly  t„;i,w  the  grand-scale  average  growth 
of  the  hemisphere  in  free  air  above  t:  e ground  surface,  as  idealized  in  Figure  4. 

While  some  lag  may  occur  due  to  drag  in  the  boundary  layer,  the  corresponding  error 
in  lead-time  is  not  nearly  as  severe  as  the  pressure  reduction  from  the  same  layer. 
Field  measurements  will  test  whether  t ,is  expectation  of  less  scatter  is  realized. 

We  also  expect  the  lead-time  will  better  "remember"  the  early  history  of  the  wave. 

If  the  explosion  starts  at  low  speed,  is  then  sustained  by  afterburning, 

the  TOA  could  be  longer,  the  lead-time  less,  than  by  an  instantaneous  explosion. 

On  the  other  hand,  compensation  occurs:  the  energetic  shock  wastes  mure  energy  early 
and  after  running  a long  time  more  slowly,  may  arrive  later  than  the  afterburner. 

UTE  form  factor  calculations  will  show  which  is  the  stronger  effect  and  how  much. 
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FIGURE  2.  LEAD  TIME  METHOD 


FIGURE  3.  LOCAL  PRESSURE  VARIATIONS  IN  SPACE  AND  TIME 
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4.  Machine  Code  and  Printout 

Appendix  A lists  the  complete  code,  is  well-annotated  and  may  be  self-explanatory. 
Written  in  advanced  BASIC  for  a personal  computer,  it  uses  two-digit  variables, 
compound  statements  in  lines  etc.,  but  in  principle  could  be  copied  and  run  as  is. 

1 nput  parameters  appear  at  the  top  of  the,  printout  (Tables  1,  2).  They  include: 

Trial  yield  Y , here  in  joules;  4.188  10  joules  = 1 KT,  4.188  10  = 1 KG, 

Mass  M of  explosive  and  immediate  case,  their  specific  energy  H relative  to  air, 
Initial  radius  R0:  HE  charge,  isothermal  sphere  (nuclear),  isobaric  sphere  (gas), 
Ambient  conditions:  pressure  P , sound  speed  C , and/or  density  Dq, 

Afterburning  fraction  AB  of  Y : Transition  pressure  P^,  strong  to  weak  shock. 

If  any  above  are  missing,  the  code  usually  supplies  a default  value  cr  computes  one. 
Input  measurements  for  evaluation  of  yield  include  pressure  P,  range  X,  and/or  TOA. 

Major  Options  are  included  regarding  input  parameters,  input  data,  and  print-out. 

Input  any  3 of  4:  yield  Y0,  initial  pressure  PI,  mass  Mi  or  specific  energy  H. 

If  any  are  omitted  (usually  P.  is  unknown),  the  code  will  calculate  it. 

If  all  are  given,  the  code  wifi  recompute  H to  make  it  consistent  with  the  rest. 

If  P,  X or  T data  are  omitted;  the  code  prints  predictions  anyway. 

Time  may  be  input  either  as  discrete  data  or  by  a fitted  curve.  (So  could  distance). 
Print-outs  can  be  predictions  only,  + diagnostics,  or  yield  analyses  with  graphs. 

Major  Computation  Blocks  (15  in  all)  are  set  off  by  remarks  in  the  LIST. 

Of  interest  as  a guide  to  the  code,  they  are  listed  on  the  first  page  of  Appendix  A. 

Computation  Procedure.  After  the  predictions  for  conditions  at  the  change  radius: 

1.  Select  a new  pressure;  next  data  or  reducing  the  previous  P by  10'*  = 1/1.26. 

2.  Compute  waste  heat  Q,  afterburning  increment  YA,  available  energy  fraction  A, 
form  factor  w/A  as  AF,  shock  velocity  U;  all  are  functions  of  the  pressure. 

3.  Compute  yield  decrement  Yl,  afterburning  YA,,  new  yield,  Y = Y.  - Y1  +YA, 
then  iterate  for  "mass-corrected"  radius  Z =(R'5  + H*Mr'  , then  get  R from  Z. 

4.  Calculate  TOA  from  dT  = dR/U,  and  new  TOA  = old  TOA  + dT.  cubed 

5.  Calculate  relative  yield  from  range,  essentially  as  (measured/calculated)cu  ® .. 

6.  Calculate  relative  yield  from  TOA  or  lead-time,  as  (measured/calculated)cu  e . 

7.  Summarize  with  an  average  yield  relative  to  input  yield,  for  range  and/or  TOA 
including  the  standard  deviation  of  the  measured  yields  about  their  mean. 

Table  1 illustrates  a printout  of  close-in  predictions  of  a massless  explosion:  Mi=0. 

The  initial  pressure  PI  and  radius  R0  mark  the  end  of  a nuclear  radiative  phase. 

The  isothermal  sphere  implies  a "square  wave",  i.e.  larger  form  factor  than  normal, 
nor  can  the  interior  gas  be  accelerated  instantly  to  reproduce  a normal  blast  wave. 
The  code  allows  it  to  do  so  gradually  by  computing  an  "inertial  mass"  as  shown. 

The  point  is:  In  a gaseous  explosion,  spark  gap,  or  any  other  non-ideal  blast  wave, 
a like  initial  dissimilarity  occurs  and  will  be  so  accommodated  with  all  UTE  codes. 

We  note  that  the  inertial  mass  found,  544  kg  = mass  of  air  engulfed  at  that  radius. 

Table  2 illustrates  another  printout  option  that  graphs  the  relative  yield  from  TOA. 

It  compares  predicted  lead-times  for  a Mk  48  torpedo  with  the  field  test  data. 

Initial  yield  included  a ground  reflection  factor  1.5  and  afterburning  of  P8XN  103. 

The  ambient  conditions  are  for  the  test  site  at  Socorro,  NM,  altitude  = 5200  feet. 

The  input  mass  was  1038  pounds,  essentially  the  warhead,  most  of  which  is  PBXN-103. 
This  was  the  very  first  test  of  the  code.  The  relative  yield  is  plotted  as  T (time). 

The  TOA  yield,  .90  predicted,  means  lead-time  itself  is  within  ,9X/  =.965,  3.5%. 

Also  shown,  not  plotted,  are  the  pressure  results:  relative  yield  1.099,  3.2%  in  R. 
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TABLE  1 
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UNIFIED  TVECRY  CF  EXPLDSICN5  (UTE),  R>M  FACTOR  NETHX>  ATO  TINES 
NJOEAR  COvfOSITE  LATA,  NXTR  72-209 


Total  Yield 

= 4. 18879E+ 12 

Input  mass  = 0 

Input  H = .25 

Ambient  pressure  = 100000 

Ambient  density 

= 1.16271 

Ambient  sound  speed  = 347 

Inertial  mass  = 

544.037 

Afterburning 

fraction  = 0 

Complete  at  Pt/Po  = 1.99526 

O’RRESSlPE 

RADIUS 

YIELD 

ARRIVAL  TINE 

LEAD  TINE 

79432.9 

4.21063 

. 839084* 

9. 50147E-06 

.0121249 

63095.8 

4.51563 

.818802 

1. 31606E-05 

.0130002 

50118.7 

4.83647 

.798519 

1 . 74959E-05 

.0139205 

39810.7 

5.17656 

.778182 

2. 26698E-05 

.0148954 

31622.7 

5.60958 

.75347 

3.00725E-05 

.0161359 

25118.8 

6.06141 

.728707 

3. 8739E-05 

.0174293 

19952.5 

6.53499 

.703973 

4. 89312E-05 

.0187839 

15848.9 

7.03314 

.679343 

6.09603E-05 

.0202075 

13600 

7.37917 

.663073 

7. 01477E-05 

.0211955 

12589.3 

7.55862 

.654897 

7. 51919E-05 

.0217076 

10000 

8.11427 

.630716 

9. 20832E-05 

.023292 

7943.29 

8.70291 

.60686 

1. 12161E-04 

.0249683 

6309.56 

9.32753 

.583393 

1. 36065E-04 

.0267445 

5011.86 

9.99126 

.560374 

1. 64565E-04 

.0286287 

3981.06 

10.6957 

.537912 

1. 98504E-04 

.0306249 

3750 

10.8641 

.532836 

2.07216E-04 

.0311015 

3162.28 

11.3621 

.518595 

2. 34559E-04 

.0325092 

2511.89 

12.079 

.499882 

2. 78187E-04 

.0345317 

1995.26 

12.8509 

.481772 

3. 30994E-04 

. 0367034 

1584.89 

13.6825 

.464255 

3.9494E-04 

.0390359 

TABLE  '2 

RUN 


UNIFIED  THEORY  CF  EXPLOSIONS  (UTE),  ROM  FACTOR  NETHX)  AND  TINES 
INHIBITOR  TESTS  ON  NK48  TORPEDCES,  JUNE  1982,  R9ELIMIINARY  DATA  SHUT  #1 
Total  Yield  = 3.85452E+09  Input  mass  = 471.252  Input  H = . 

Ambient  pressure  = 83000  Ambient  density  = .991783 

Ambient  sound  speed  = 342.29  Calculated  Initial  Pressure  = 1.08279E+07 


Afterburning  fraction  = .3  Complete  at  Pt/Po  = 1.99526 
P/PO  TQ4  Yield  Nfeas.  T .6  Relative  TOA  Yield 


.456874 

. 884941 

.0679804  . 

T 

.376298 

. 762048 

.0809575  . 

T 

.223453 

.938644 

.128332  . 

T 

.201855 

.88243 

.140158  . 

T 

.154507 

.983428 

.176062  . 

T 

.0930362 

.933178 

.268408  . 

T 

.0880522 

.915598 

.280596  . 

T 

1.4 


Yield,  relative  to  input  = 1.09959 

Standard  deviation,  % = 21.1099  based  on  7 samples,  P=>  1 psi 
TOA  Yield,  relative  to  input  = .900058 

Standard  deviation,  TOA  yield,  % =•  7.77284  based  on  7 samples 
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v**r''* 


5.  Test  of  Methods  with  Nuclear  and  HE  Data 

The  form  factor  and  lead-time  methods  were  tested  against  a broad  spectrum  of  data. 
Nuclear  data  check  on  absolute  yields  by  their  radiochemical  and  hydrodynamic  yields 
and  check  the  equation  of  state  more  severely  at  higher  pressures  than  HE  reaches. 

HE  data  check  non-ideal  effects  like  large  mass,  afterburning,  and  secondary  shocks. 
The  broad  range  of  data  checks  for  self-consistency  and  exposes  systematic  errors. 
Figures  6 to  iO  graph  the  detailed  results  and  Table  3 summarizes  them. 

Blast  theory  is  usually  checked  against  data  by  pressure-distance  plots  like  Fig.  3. 

But  as  seen  there,  UTE  matches  composite  data  so  closely  that  graphs  are  inadequate. 
Instead,  UTE  computes  the  relative  yield  at  each  pressure  level  and  we  plot  that. 

On  Figure  5,  the  line  widths  approach  3%  in  radius,  10%  in  yield,  too  small  to  see. 

On  Figures  6 to  10,  the  three  central  lines  are  relative  yields  of  1.0  - 10%, 

as  if  the  graph  of  Fig  5 were  blown  up  to  broaden  the  lines  to  the  band  width  shown. 

17  5 

1 K T Nuclear  Composite  ’ (Fig  6)  covers  from  13600  to  .07  bars,  10  times.  ^ 

The  average  yield  1.024  KT  - 5.3%  matches  the  line  width  of  the  source  curve  ^ 

and  is  significant  because  the  high  pressures  are  superbly  accurate  fireball  data. 

The  TOA  yields  also  are  excellent  at  high  pressures;  overall  is  1.08  KT  - 14%. 

The  excursion  at  low  pressures  is  probably  due  to  a fitted  time-of -arrival  curve. 

The  consistency  in  yield  is  assurance  that  the  high  pressure  UTE  equation  of  state 

is  realistic  relative  to  the  ideal  gas  law,  used  for  air  at  pressures  below  10  bars. 

KING  Fireball'*  (Fig  7)  is  probably  the  best  pressure-distance  data  in  existence: 
high  yield,  air  drop,  negligible  mass  effect,  all-fission,  a perfect  circle  fireball. 
Radiochemistry  gave  545  KT,  hydrodynamic  yield  555  KT,  fireball  scaling  595  KT. 

Here,  pressure  and  TOA  both  give  586  KT;  scatter  of  3-7%  is  round-off  error  in  data. 
This  one  test  is  definitve:  all  the  measured  data  are  digital  —no  graphing  errors. 

O 

Nuclear  Blast  Standard  (Fig  8)  is  not  data  but  an  artifical  viscosity  hydrocode. 

The  absolute  value  of  yield  .997  KT  checks  superbly,  but  the  scatter  is  over  14%. 

Its  initial  pressures  are  known  to  be  50%  lowborn  actual  fireball  theory ^a^  data. 

At  low  pressure  its  P-R  curve  decays  like  R ’ , flatter  than  UTE,  P~R  * . 

Classic  theory  predicts  R"  , but  field  measurements  always  decay  much  faster. 

Q 

1 KG  TNT  Composite  (Fig  9):  splendid  agreement/consistency,  up  to  the  charg^,^ 
and  for  .07<P/P  <2,  the  UTE  calculation  agrees  well  with  often-measured  P~R 
The  excursion  below  .07  bars  is  probably  not  real,  but  old  data  piously  fit  to  R . 

The  absolute  yield  is  .714  KG,  714  cal/gm;  earlier  UTE  methods  gave  720  cal/gm. 

o 

1KG  H6  Composite  (Fig  10  ) is  a check  with  a heavily  aluiminized  explosive, 
where  the  afterburning  fraction  is  estimated  as  .30.  The  consistency  5%  is  suberb. 

The  absolute  yield  is  1.014  KG  HE,  or  1014  cal/g,  equivalent  to  1.4  times  TNT. 

Previous  UTE,  DSC*  (not  shown)  has  been  used  successfully  on  so  many  NESIP  and 
other  tests  that  it  is  of  interest  to  use  a DSC  calculation  (M=0)  as  input  data  here. 
The  result:  Relative  yield  1.00000,  - 3.5%,  no  sensible  difference  between  them. 
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FIGURE  6.  RELATIVE  YIELD  VS  OVERPRESSURE,  1 KT  NUCLEAR  COMPOSITE 
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OVERPRESSURE  RATIO  P/PD 
FIGURE  10.  RELATIVE  YIELD  VS  OVERPRESSURE,  H-6 


6.  MK.  48  Torpedo  Inhibitor  Tests 

2 10 

Recent  field  experiments  on  the  design  of  inhibitors  ’ for  the  Mk  48  torpedo 
provide  an  opportunity  to  test/apply  these  new  methods  to  a typical  NES1P  problem. 

The  test  set-up  is  shown  at  the  top  of  Table  4j  essentially,  it  was  a donor  warhead 
flanked  by  two  acceptor  warheads,  with  inhibitor  plates  between  them,  on  each  side. 

On  one  side  the  plates  were  steel,  on  the  other  side  aluminum,  thicknesses  as  shown. 
Detonation  was  suitably  instrumented  with  witness  plates,  cameras  and  flash  panels 
and  in  a blast  line  by  pairs  of  pressure-time  gages  near  140,  240,  3U0,  400  feet. 

Even  so,  critical  questions  arise  in  all  such  tests: 

If  the  witness  plate  did  hole,  did  full  detonation  of  the  entire  warhead  occur? 

If  it  did  not  hole,  could  the  charge  have  moved  and  a delayed  detonation  occur? 

To  understand  either  case,  we  need  to  measure  the  yield  output  for  each  event. 

The  test  results  were  unequivocal  and  corroborative  among  all  the  test  evidence. 

Shots  1 and  2:  no  acceptor  detonated.  Shot  3:  the  aluminum  side  holed  in  situ. 

The  140  and  240  foot  gages  on  Shot  //3  had  double  pulses  that  coalesced  by  300  feet. 
Still  the  critical  questions  remain:  How  much  energy  did  each  shot  yield? 

Peak  pressure  results  are  compared  on  Figure  11  with  the  pre-shot  calculations. 

The  data  on  //I  seem  somewhat  low,  on  //2  somewhat  high.  But  scatter  makes  it  doubtful: 
Excepting  t.wo  "low"  points  on  #1,  one  "high"  on  //2,  13  remaining  points  replicate. 

Shot  //3  leaves  no  doubt  the  curve  beyond  300  feet  represents  twice  the  yield. 

The  corresponding  pre-shot  estimates  and  lead-time  data  are  shown  on  figure  12. 

Now  there  is  no  doubt  that  Shot  // 2 was  larger  than  #1,  nor  that  Shot  //3  was  double. 
Considering  this  was  the  first  test  of  a lead-time  prediction  on  HE,  it  checks  well. 

Relative  yields  on  shots  ) and  2 are  plotted  on  Figures  13-14,  summarized  on  Table  4. 
Compare  the  pressure  results:  1012  KG  - 21.1%  vs.  1386  KG  - 17.9%. 

The  ratio  1386/1012  - 1.37  is  impressive,  except  that  37%  is  not  far  different  from 
the  arithmetic  or  the  Pythagorean  sum  of  deviations,  21.1+17.9.  One  is  just  not  sure. 

Now  compare  lead-time  yields:  828  KG  - 7.8%  vs.  1108  KG  - 6.B%. 

Again:  1108/828  = 1.34  is  impressive  and  34  is  more  than  twice  any  sum  of  7.6  and  6.8 
These  confidence  levels  make  a strong,  objective  case  for  the  merits  of  lead-time. 

As  Table  4 indicates,  the  yields  on  shot  3 were  definitely  doubled,  by  either  method, 
but  no  predictions  had  been  made  with  history  effects  for  catch-up  of  second  pulses. 

We  have  yet  to  resolve  why  the  lead-time  gives  lower  yields  on  both  shots  1 and  2. 
Compare  range/lead-time:  1012/828  = 1.23  (shot  1)  and  1306/1106=  1.25  (shot  2). 

It  is  precarious  to  prognosticate  with  preliminary  data  until  they  really  do  firm  up, 
but  two  main  ideas  are  noteworthy  here:  1.  measuring  sound  velocity,  2.  reflection  . 
Among  many  ways  to  measure  sound  velocity  —absolute  temperature  + wind  velocity, 
a microcharge  fired  just  before  the  main  shock,  or  compute  C from  the  P-t  data— 
all  three  differed  at  the  field  tests,  and  we  have  not  yet  resolved  why. 

On  the  other  hand,  the  lead-time  could  well  be  telling  us  a real  fact: 

The  calculations  assumed  the  torpedo  explosion  reflected  instantly  off  the  ground. 

Whereas  we  know  it  must  have  run  for  some  time  as  a free  air  explosion, 
and  was  slowed  by  the  inhibitor  plates  and  by  the  acceptor  in  that  direction. 

We  also  know  that  the  shock  is  slowed  by  the  dust-load  in  the  boundary  layer 
The  present  results  are  based  on  sound  velocity  as  calculated  from  pressure  gages. 
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TABLE  4. 

MK  46  TORPEDO  INHIBITOR  TESTS 


DONOR 


YIELD  (REFLE  1.)  CONSISTENCY 


SHOT  =1  1 1/2"  ALUMINUM 

NOMINAL 
1000  KG 

SPAN,  PSI 
10  * 1 

RANGE:  1012  KG 

' 21.1% 

=2  3/4”  ALUMINUM 

1000  KG 

10  -*  1 

TOA:  828  KG 
RANGE:  1386  KG 

* 7.8% 

• 17.9% 

=3  3/8"  ALUMINUM 

1000  KG 

10  • 1 

TOA:.  1108  KG 
DOUBLED. 

* 6.8% 

DOUBLED: 

PRELIMINARY  DATA  AND  ANALYSES 


r 
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LEAD  TIME,  SECONDS 


FIGURE  1Z  TEST  OF  LEAD  TIME  OF  MK  48  TESTS 
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MK  48  TORPEDO  INHIBITOR  TESTS 
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7.  Conclusions 


1.  Form  Factor.  The  form  factor  method  with  the  unified  theory  of  explosions  (UTE) 
agrees  well  with  nuclear  and  high  explosive  data,  with  earlier  methods  of  UTE 

and  offers  a facile  way  to  describe  non-ideal  and  non-spherical  explosions. 

2.  Lead-Time.  The  lead-time  method  is  a simply  instrumented  way  to  measure  yield 
at  high  and  low  shock  strengths,  with  much  less  scatter  than  pressure  measurements. 

3.  UTEFORM.  Form  factor  and  lead-time  together  offer  a new  powerful  diagnostic  tool 
to  solve  the  unpredictably  broad  problems  which  explosion  safety  requires  such  as 
sympathetic  reactions,  early  blast  history,  unusual  afterburning  or  energy  release. 

4.  Absolute  Yield.  The  definition  10 ^ cal/KT  = 106  cal/KG  = 10^  cal/gm 

is  a modern  rational  way  to  correlate  any  explosion:  nuclear,  chemical,  other  source. 

It  is  necessary  because: 

Different  hE's  do  not  necessarily  scale  with  each  other  nor  with  other  sources. 

That  is,  equivalent  weight  is  certainly  not  constant  at  high  shock  strengths 
and  is  not  necessarily  constant  even  at  acoustic  shock  strengths. 
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Appendix  A 


A 


LIST  for  UTEFORM 

Form  Factor  and  Lead-Time  Methods  With  the  Unified  Theory  of  Explosions 

Major  Computation  Blocks  by  Line  Numbers: 

0-199  Input  parameters 

200-299  Compute  PI.  given  Y0  and  H 

299-300  Compute  M or  H,  given  Y0,  PI  and  either  M or  H 

400-450  Compute  trial  Y,  if  unspecified 

450-500  Print  column  headings 

500-570  Data  processing  and  pressure  selection 

570-699  Energy  gains,  losses,  new  yield  and  range 

700-800  Equation  of  state  sub-routine 

800-850  Energy  gain  and  loss  sub-routines 

850-899  Form  factor  sub-routines 

900-999  Time-of-Arrival  sub-rountine 

1000-1400  Input  data:  pressure,  distance,  time 

1400-1500  Example:  Fitted  time-of -arrival  sub-routine  for  a nuclear  composite. 
1500-1599  Relative  yields  from  range  and  time-of -arrival 
1600-2000  Yield,  standard  deviations  and  termination. 


Hints: 

1.  Any  consistent  set  of  units  may  be  used.  ^ 

If  energy  is  in  joules,  R in  meters,  M in  kg,  then  is  in  pascals  =10”  bars. 
Line  50,  as  written,  converts  from  KT  (4  pi/3  *10  ; to  joules; 

use  line  51  to  enter  the  KT,  KG  or  cal/gm. 

2.  Chanqe  data  with  a line  editor,  it  will  save  retyping  the  remarks  in  that  line. 

3.  In  general,  the  variables  are  defined  by  remarks  the  first  time  used  in  the  LIST. 

4.  For  help,  call  Fran  Porzel,  202  394  1166  (office)  or  703  533  7973  (home). 


rj  w ii*  in Wl - i ’ 


0 J =2  '0  = Predictions,  1=  analysis,  2,3  = range  or  TEA  yield 

1 REM.  .Any  self  consistent  units  may  be  used;  P is  in  pascals  for  Y kg,  R meters 
5 PRINT  "M£DvP/UT8",  OATE$,  TIMES:  IF  J=0  THEN  15 

15  PRINT  "UNIFIED  THECRY  CF  EXPLOSIONS  (LfTE),  FCFM  FACTOR  NCTVQD  AO  TItvES" 

23  PRINT  "NUCLEAR  COMPOSITE  LATA,  NXTR  72-209" 

25  B = 4*3.141592/3  'Form  factor  for  sphere;  mech.  eq  of  heat,  13.6  deg  cal. 

30  P0  = 1E5  'Ambient  pressure;  1 bar  = 1E5  pascals  =1E5  kg/m/sec*2 

31  P0  = 1*P0  'Erases  possible  previous  entry 

35  £0  = 2.5  : K0  = 1/EG  +1  'energy  and  adiabatic  coefficients,  ambient 

40  DO  =1.1613  :C0  = 1138. 45*. 3048  'input  ambient  density  DO  or  sound  speed  CO  or  both 

41  DO  = 1.11613  'Erases  possible  previous  entry  for  Do 

43  IF  CO  > 0 THEN  DO  = K0*PQ/CG*2  : CDTO  50  'Override  DO  by  equation  of  state 
46  IF  CO  = 0 TVEN  CO  = SGR(K0*P0/DO) 


50  Y0  = 8*10*12 

51  Y0  =Y0*1.0 

52 

53  ABr.OO 

56  Y0  = Y0*(1 -AB) 

60  R0  =4.2 

61  R0  = 1*R0 
70  H = .25 
80  M =0 

83  M = H*M0/B/D0 

85  IF  RO  = 0 TVEN  RQ 

86  ZO 
90  PI  = 8E9 


'Yield;  1 KT  = 10*12  cal  = 4pi/3*lel2  kg  m*2/m*3/sec*2 
'Relative  yields  from  earlier  runs  or  fits 
' 1 KG  = 10*6  cal  = 4pi/3*le6  kg  m*2/m*3/sec*2 

'Afterburning  fraction 
'Yield  before  afterburning 

'Radius  of  isothermal  sphere  or  charge  radius  , 1 KT 
'Erases  possible  previous  entry 
'specific  energy  of  debris  to  air;  use  .5  for  massive 

'converts  mass  to  equivalent  voltme  of  air 
(M3/B/1500)*(l/3)  ELSE  86  'Replace  1500  w/  D of  charge 
(R0*3  + M)*(l/3)  'Z  = Sshock  radius  corrected  for  UTE  mass  effect 

93  PT  = PG*10* . 3 'Transition  pressure,  bock-keep  end  of  afterburning 

95  QZ  = 3.5:  YZ  = .5  'Default  intial  values  for  dlnQ/dlnZ,  dlnY/dlnZ  (ideal) 
100  PRINT  "Total  Yield  ="Y0/(1-AB),  "Input  mass  ='TvO,  "Input  H ="H 
120  PRINT  'Ambient  pressure  r"P0,  "Ambient  density  ="D0 
140  PRINT  "Ambient  sound  speed  ="C0, 

199  REM CPTION  TO  CALCULATE  PI , GIVEN  Y0  AND  H 

200  P = PI  : IF  PI  >0  TVEN  300 

210  PI  = YO/B/ZQ'3  'Trial  value;  A*F  approx  1 for  strong  shocks 
220  P = PI 
230  CDSIB  700 
240  AF  = A/3 
250  Y = Y0*( 1-Q/P) 

260  P = Y/B/AF/Z0*3 
270  IF  ABSCP/PI  -1)  <,001  TVEN  285 
280  PI  = P:  GOTO  230 
285  PC  = PI 


'Estimate  for  waste  heat  of  charge  or  isothermal  sphere 


Iterate  for  PI 

'Save  revised  pressure  PC  at  charge  surface 


'Waste  heat  in  radiative  phase  or  chargew/af  = 1 


290  PRINT  'Calculated  Initial  Pressure  ="P 
296  GOTO  445 

299  REM CPTICN  TO  CALCULATE  H,  GIVEN  PI  AND  Y. 

300  CDSUB  700 
315  IF  Y0  =0  TVEN  400 
320  Y = Y0*(  1-Q/P) 

330  IF  M3>0  AND  H>0  TVEN  390 
340  ZQ=(Y/B/PI/AF)*(l/3) 

350  MH  = 8*D0*(Z0*3  -R0‘3) 

360  IF  M0  = 0 TVEN  M3  = NH/H  ELSE  380 

370  PRINT  "Inertial  mass  ='TvC  : IF  3 = 0 THEN  380 

380  IF  H = 0 TVEN  H = NH/MO  ELSE  390 

385  RRINT  "Calculated  specific  energy  H ="H 

390  M = VHvC/8/D0  'Computes  "inertial  volume"  from  apparent  mass 


330 


399  REM OPTION  TO  CALCULATE  Y (Needs  debugging) 

400  IF  Y = 0 THEN  Y = B*PI*A*F*Z0'3  ELSE  430 
405  GOTO  930 

410  Y0  = Y/(l-  PI/P0r(l/K  -1)) 

420  PRINT  "Calculated  initial  yield  YQ  ="Y0 
430  P = PI 

445  PRINT"Afterburning  fraction  ="AB,  "Complete  at  Pt/Po  ="PT/P0 

450  IF  J<2  THEN  PRINT"0'PRESSURE",  "RADIUS",  "YIELD", "ARRIVAL  TIME", "LEAD  TIME" 

470  IF  J<2  THEN  PRINT  " ",  "Measured", "Range  Yield", "Measured",  "TOA  Yield" 

475  IF  J = 1 THEN  PRINT  " Q/P","  Z",  " YZ",  " QZ".  " AF" 

480  IF  J = 2 THEN  PRINT  "O'Press."TAB(10)"Rel.  Yld"  TAB(20)"dlnY/dlnZ"  TAB(30)".6", 

490  IF  J = 3 THEN  PRINT  "P/P0"  TAB(10)"TOA  Yield"  TA6(20)"Meas.  T"  TAB(30)".6", 

491  F J = 2 OR  3 THEN  PRINT  "Relative  Yield"  TAB(69)"1.4" 

499  REM PRESSURE  SELECTION 

500  READ  PX,X  : IF  PX  = 0 THEN  1600 
510  PX=  PX*1£5 

513  REM  Use  GO'i 525  with  no  measured  times  or  if  TX  is  in  seconds 
515  GOSU0  1400  'Sub-routine  for  fiited  TOA  curve 

520  TX  = TX/1000  'Fitted  curve  was  in  miitiseconds 

525  IF  P>PX  THEN  N = .1*1NT(10*LOG(PI)/LOG(10)  ) : GOTO  540 

526  PI  = PX  H=0  ; X0  =X  : GOTO  300 
530  N = N - 0.1 

540  P = 10 "N 

545  IF  ABSCPI/P  -1X.001  THEN  530 
550  IF  P <PX  THEN  P = PX 
560  GOSUB  700 

570  REM, GAINS,  LOSSFS,  NEW  YIELD  AND  ARNGE 

580  IF  QI  = 0 THEN  630  'By-pass  gain  or  loss  at  RO;  avoid  /Q  error  @600 

590  YH  = Y +YA  'Add  afterburning;  hold  Y as  Z is  iterated 

595  F QZ  = 0 THEN  605  'IF  Q = QI  THEN  Z =ZI  and  Y1  = 0 
600  Y1  = 3*B*QI*(Zr3)/(aZ  -3)*(1  - (Q/QIXU-3/QZ)) 

605  IF  P<PT  THEN  YA  = 0 'No  afterburning  beyond  transition  pressure 

610  Y = YH  -Y1 

630  Z =(Y/B/P/AF)*(l/3) 

631  IF  QI  = 0 THEN  637 

632  IF  Q = QI  OR  Z = ZI  THEN  635  'Circumvents  repeateu  pressure  problem 

633  QZ  = LOG(QI/Q)/LOG(Z/ZI) 

635  IF  ABS(Z2/Z  - DX0U001  THEN  Z2  = Z : GOTO  595 
637  IF  M>Z‘3  THEN  Z = (M  *-R0‘3rU/3) 

640  R = (Z~3  - MFd/3) 

645  ON  ERROR  GOTO  650  'Avoids  d/0  on  initial  pass 

646  YZ  = LOG(Yl/Y)/LOG(Z/ZI)  'Calculate  dlnY/dlnZ  for  later  use  in  870 

650  GOSUB  900.  ’Get  time  increment 

660  T = T +T1 
675  IT  J >1  THEN  687 
680  PRINT  P/P0,  R,  Y/Y0,  T,  R/C0  -T 
687  IF  ABS(P/PX-1)<.001  THFN  GOSUB  1500 
690  IF  J = 1 THEN  PRINT"  "Q/P,  Z,  " "YZ,  " "QZ,  " "AF 

696  PI=P  : QI  = Q : ZI  =Z  : RI  =R  rUI  =US  :YI  = Y:AI=AF 
698  IF  ABS(P/PX  -1X.001  THEN  500  ELSE  530 
Ready 
> 
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699  REM EQUATION  CF  STATE  SUB-FOLfTINE. 

700  P = P/PO 
710  IF  P<10 
720  IF  PC100 
730  IF  P<700 
740  IF  PC1000 
743  IF  P<4000 
750  IF  P<40000 
755  IF  P540000 


'Equation  of  state  is  described  best  by  pressure  ratio 
THIN  £=2.5  :CDTO  760  'E  is  ratio  =energy/PV,  i.e.epsilonh 
THEN  E=  2.5  + 1. 5*LDG(P/10)/LQG(  10)  :CDTO  760 

OEN  E=  4.0  + 1.55*LGG(P/100)/LOG(7)  :ODTO  760 
THEM  E =5.55  -. 55*LGG(P/ 700) /LOG(  10/7)  :CDTO  760 
THEM  E =5.0  - LCG(P/1000)/LCG(4)  :CDTD  n6Q 

OEN  E = 4.0  :GDTO  760 

OEN  E = 4.0  - . 67*LCG(P/40000 ) /LQG(  2 ) :GOTO  760 
760  O = ((2*L  +1)*(1  +P)  +1)/(P  + 2*£0  +2)  ' density  ratio,  real  gas 

770  A = E*(  1 +P)/P*(1-  (1+P)"(-1/(E  +1)))  'Pronrpt  energy  factor  A =Vifork/(P-P0)V 
780  K = 1/E  +1  'Same  as  epsilon  = l/(k-l) 

790  P = P*P0  'get  back  to  absolute  overpressures 

799  REM ENERGY"  LOSS  AND  GAIN  SUBROUTINE 

800  IF  P/P0  < .06  THEN  840  'Q  wi  1 1 soon  truncate  to  0 if  you  don't  do  this 

805  IF  P/P0  <11.3  TVEN  830  'Exact  match  w/  ideal  gas  i P=11.3,  3.4 

810  L = .4342948*LOG(P/P0)  'convert  pressure  ratio  to  log  base  10 

820  Q = P0*HT((21.75-L)*(L-1)/16)  :GDTO  845  ’Semi -empi r ica 1 fit  for  reel  air 
830  Q = P0*((1  +P/P0r(l/K)/D  - 1 ) / (K-l ) :CDTO  842  'classical  adiabat 

840  Q =P0*(K+1 )*( (P/P0/K) * 3 )*(  1 - 1.5*P/P0)/12  'acoustic  dissipation 

842  IF  P>PT  TVEN  845  'Argument:  wave  form  and  losses  are  manifest  at  shock 

843  Q = Q*(l  +AB)  'Argument:  secondary  shock,  other  losses 

845  ZP=  1/4  'dlnZ/dlnP;  assunes  VA  goes  as  volume  and  timt 

846  YA  = AB*Y0*((PG/PrZP  -(PC/PI )‘ZP)/((RE/PT)*2P  - 1)  'A8  is  prop,  to  Z-zi 

847  IF  P<PT  OEM  YA  = 0 

850  REM AF  = A*F  SUBROUTINE 

860  IF  P>PT  OEM  AF  =A/3/(l  +P0/P)*2/(1  + A8)  :GDTO  890 

861  REM:  Strong  shock,  F=  .42,  mean  A=. 8*A( shock ) , /(I  +AB)  is  peaked  wave  form 
865  IF  YZ  = 0 OEN  YZ  = 1 

870  IF  P<PT  THEM  AF  =3*Q/P/YZ  'Weak  and  second  shock,  YZ  stable 

880  IF  AF> AI  CR  AF<=0  OEN  AF=  AI  'By  passes  troubles  at  transition  pressure 

890  RETIRN 

900  REM TINE-CF  ARRIVAL  SUB-ROUTINE 

910  IF  P/P0  >450  OEM  K =1.2  + .2*LGG(P/P0/450)/LGG(  2)  : GOTO  940  'real  gas  K 
920  IF  P/P0  >15  OEN  K = 1.4  - .2*U0G(P/P0/15)/LOG(30):  CDTO  940  'real  gas  K 

930  K = K0 

940  US  = SQR(P/D0/(1  -1/D))  'used  previous  K to  calculate  U 

945  IF  Zl  = 0 OEN  UI  = US  ' UI  not  yet  initialized  as  in  696 

950  UB  = (1/US  + l/UD/2  'Mean  for  integrating  dt  as  dX/U 

960  T1  = UB*(R  -RI)  'Time  increment 


970  IF  ZI  = 

0 AND  P/P0 

>10000  OEM  T1 

= .2*T1  : GOTO  990  'Radiative  phase 

980  IF  ZI  = 
990  RETURN 

0 OEN  T1  - 

: . 5*T1 

'ball  park  estimate  for  detonations 

999  lEM. . . . 

1000  CATA  13600,  7,32, 

3750,  10.7, 

1550,  13.7 

1010  DATA 

1000,  15.5, 

510,  19.2, 

200,  25.6 

1020  DATA 

100,  32.2, 

50,  41, 

20,  57.3 

1030  DATA 

10.,  75.1, 

8,  82.4, 

6,  91.5 

1040  DATA 

5,  98.8, 

4,  108.5, 

3,  122.5 

1050  CATA 

2,  147.4, 

1.0,  208, 

.5,  302 

1060  DATA 

.2,  544, 

.10,  905, 

.07,  1200 

1400  REM SLB -ROUTINE  FCR  FITTED  TINE  CF  ARRIVAL 

1410  IF  (PX/PO)  < 500  THEN  1430 

1420  TX  = (PX/PO/ 500 ) ' ( -4/ 5 ) r QOTD  1490 

1430  IX  = (PX/PO/ 500)' (-6/7) 

1435  TX  = TX*TX'(LOG(TX)/175 

1440  IF  A8S(T2/TX  -1)  < .0000.  THEN  1490 

1450  T2  = TX  :CDTO  1430 

1490  RETURN 

1499  ITEM RANX  AND  TINE  CF  ARRIVAL  YIELD 

1500  IF  A8S(P/PX  -1)>,001  THEN  1595  'Passes  only  measured  points 
1510  YX=  (XA3  hM)/2'3 

1533  TB  = INT(50*YX  + .25) 

1534  IF  TB>77  THEN  TB  = 77 

1535  IF  J = 2 THEN  PRINT  P/PO  TAB{10 )YX  TAB(20)YZ  TAB( 30)".  ”TA8(TB)"Y"  TAB(70)"." 

1545  IF  US /CO  >2  THEN  YT  = (TX/T)A3  •.  CD  TO  1555 

1550  YT  = ((X/CO  -TX)/(R/C0  -T))*3 
1555  TT  = INT ( 50*YT  + .25) 

1560  IF  TB>77  THEN  TB=  77 

1570  IF  J = 3 THEN  PRINT  P/PO  TAB(  10 )YT  TAB(20)T  TAB(30)"."  TAB(TT)"T"  TAB(70)"." 

1575  IF  P<.068*P0  THEN  1590 

1576  1=1+1 
1580  SX  = SX  +YX 

1585  VX  = VX  + YX'2 

1586  IT  = IT  +1 

1587  ST  = ST  + YT 

1588  VT  = VT  + YT‘2 

1590  IF  J = 1 THEN  PRINT"  Meas:",  X,  "*"YX,  TX,  »*"YT 
1595  RETURN 

1599  REM YIELDS  AND  STANDARD  CEVIATICNS 

1600  IF  I <2  IVEN  2000 

1602  S = VX/O-l)-  (SX*2)/I/(  1-1) 

1604  IF  IT  < 2 IHEN  1610 

1606  S2  = VT/OT-l)  -(ST*2)/IT/(IT-1) 

1608  PRINT 

1610  PRINT  "Yield,  relative  to  input  ="  SX/ I 

1630  HRINT  "Standard  deviation,  % ="  100*SGR(S)*(  I/SX)  "based  on"  I "samples,  P=>  1 psi" 
1640  IF  ST  = 0 THEN  1680 

1650  PRINT  "704  Yield,  relative  to  input  ="  ST/IT 

1670  PRINT  "Standard  deviation,  104  yield,  % ="100*SGR(S2)*( IT/ST)"based  on  "IT"samples" 
1990  Q4TA  0,  0,  0 
2000  END 


